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Abstract 

One way of obtaining a version of quantum mechanics without observers, and 
thus of solving the paradoxes of quantum mechanics, is to modify the Schrodinger 
evolution by implementing spontaneous collapses of the wave function. An explicit 
model of this kind was proposed in 1986 by Ghirardi, Rimini, and Weber (GRW), 
involving a nonlinear, stochastic evolution of the wave function. We point out 
how, by focussing on the essential mathematical structure of the GRW model and 
a clear ontology, it can be generalized to (regularized) quantum field theories in a 
simple and natural way. 

PACS numbers: 03.65.Ta; 03.70.+k. Key words: quantum field theory without 
observers; Ghirardi-Rimini- Weber model; identical particles; second quantization. 

1 Introduction 

John S. Bell concluded from the quantum measurement problem that "either the wave 
function, as given by the Schrodinger equation, is not everything or it is not right" 
[I]. Let us assume, for the purpose of this paper, the second option of the alternative: 
that the Schrodinger equation should be modified in such a way that superpositions of 
macroscopically different states, as exemplified by Schrodinger's cat, either cannot arise 
or cannot persist for more than a fraction of a second. Theories of this kind have come 
to be known under the names of "dynamical reduction" , "spontaneous localization" , or 
"spontaneous wave function collapse," and have been advocated and studied by vari- 
ous authors 122 123 D33 EH 13 QUI 123 123 E 123 EE3; see for an overview. The 
merit of such theories, which I shall call "collapse theories" in the following, is that they 
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provide, instead of statements about what observers would see if they were to make 
certain experiments, a possible story about what events objectively occur: they are, in 
other words, quantum theories without observers. In collapse theories, observations are 
merely special cases of the objective events, for which the theory, if it is to be empir- 
ically adequate, predicts the same distribution of outcomes as the standard quantum 
formalism. An example of a no-collapse quantum theory without observers is Bohmian 
mechanics [TJJ |3] . 

An explicit collapse theory has been proposed by Ghirardi, Rimini, and Weber 
(GRW) in 1986 [TS] for nonrelativistic quantum mechanics of N distinguishable par- 
ticles. While the testable predictions of the GRW model differ in principle from the 
quantum formulas, as yet no experiment could decide between the two, provided one 
chooses the collapse rates proportional to the masses, since the differences are too tiny 
for standard experiments and the experiments leading to noticeable differences are hard 
to carry out In the GRW model, the wave function obeys the unitary Schrodinger 
evolution until, at an unforeseeable, random time, it changes discontinuously in an un- 
foreseeable, random way — it collapses. In order to obtain similar models for quantum 
field theories (QFTs), one path of research, which has been followed under the name 
"continuous spontaneous localization" (CSL) [2UH7J, * s based on the idea of incessant 
mild collapses, so that the quantum state vector follows a diffusion process in Hilbert 
space. For a CSL model for identical particles, see [2*lj . 

I will show how the GRW model can be extended to QFT in a more direct way; the 
resulting collapse QFT is as hard to distinguish experimentally from standard QFT as 
the GRW model from standard quantum mechanics, while solving the quantum mea- 
surement problem in the same way as the GRW model. My proposal retains (and indeed 
is based on) the discreteness of the GRW model; it gets along with indistinguishable par- 
ticles, both fermions and bosons, and with particle creation and annihilation. However, 

1 will not try here to make the theory Lorentz-invariant; I hope to be able in a future 
work to combine the construction I present here with the Lorentz-invariant version of 
the GRW model for N particles that I have described recently |2~8] . When applied to a 
system of N identical particles, my proposal yields a collapse process proposed already 
in 1995 by Dove and Squires jTTJEl], which, however, seems to have received little or 
no attention so far. 

2 Mathematical Framework of the GRW Model 

I now describe the GRW model in an unusual way that emphasizes the abstract mathe- 
matical structure it is based on and uses an ontology proposed by Bell jH E] . According 
to this flash ontology, matter consists of millions of so-called flashes, physical events 
that are mathematically represented by space-time points. The flashes can be thought 
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of as replacing the continuous particle trajectories in space-time postulated by classical 
mechanics. The flashes are random, thus forming what probabilists would call a point 
process in space-time, with a distribution determined by the (initial) wave function. 
The reader may wonder why, in a collapse theory, there is any need at all to introduce 
space-time objects such as flashes, a question that I will take up in Section For 
now I ask for patience and suggest to regard it as the sole role of the wave function to 
determine the distribution of the flashes. 

We begin with an (arbitrary) Hilbert space ffl with scalar product (<p\ip), which 
in the GRW model is L 2 (R 3N ), and an (arbitrary) self-adjoint operator H on Jif, the 
Hamiltonian, which in the GRW model is the usual Schrodinger operator — y-A + V. 
The third mathematical object that we need, besides Jif and H, is a 

positive-operator- valued function A(ac) (1) 

on physical space M 3 acting on Jif. (For mathematicians I should add the postulate that 
there is a dense domain in J$f on which all of the A(x) are defined.) A forms the link 
between Hilbert space and physical 3-space and can be thought of as representing the 
(smeared-out) position observable, and in particular as representing a preferred basis. 
In the GRW model, A (a;) is a multiplication operator 

A(aj) ...,r N )= Af e ^ x - r ^/ a2 if){r u ...,r N ), (2) 

multiplying by a Gaussian with width a/\/2 (in Bell's |I] notation) and center x, where 
Af and a are new constants of nature with suggested values of 

N w 10 5 sec _1 m" 3 , a^l(T 7 m. (3) 

We are considering for the moment only flashes correspoding to the z-th particle coor- 
dinate r*j, i G {1, . . . , N}, and will later consider flashes due to all particles. 

The role of A is to define the rate of a flash to occur, i.e., the probability per unit 
time. That is why I will call A(x) the flash rate operators. The rate at time t = of a 
flash in the set B C R 3 is 

(V|A(BM, (4) 

where A (B) is short for 

A(B) = [ d 3 a;A(a;). (5) 
Jb 

This shows why A(B) must have dimension 1/time. In the GRW model, the total flash 
rate is independent of the quantum state if) since A(M 3 ) = 7i 3 ^ 2 J\fa 3 I is a multiple of the 
identity operator I on Jif. The constant in front of J, or, equivalently, the total flash 
rate, is called 1/r in Bell's jl] notation, with 

r^l0 15 sec. (6) 
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It has been pointed out by Pearle and Squires [26] that a fixed r would contradict the 
observed stability of nucleons; they suggest instead that the flash rate constant A/a 3 be 
proportional to the mass of a particle, rather than a universal constant. In the QFT 
model we will devise, however, A(IR 3 ) will not be a multiple of the identity, and this 
is the only aspect in which we essentially generalize the mathematical structure of the 
GRW model. 

We now define the probability distribution of the first flash (Xi,Ti) with random 
location Xi and random time Ti, a probability distribution on the space-time region 
with t > 0. The distribution is quadratic in ip. The wish that the rate be given by A and 
that the evolution before the flash be given essentially by H leads us to the following 
form for the distribution: 

P(Xi G d 3 a;, T x G dt) = W? A(x)W t \ip) d 3 x dt , (7) 

where the star denotes the adjoint operator, and 

W t = e -5 A ( K3 )*-i^ for t > , W t = for t < . (8) 

Without the A(M 3 ) term, this would be the ordinary unitary evolution; we need the 
additional term to keep track of the probability that time t is reached without a flash. 
Indeed, the definition (jHJ) implies that (JJJ) is a probability distribution: 



d 3 x / dt(ij\W t *A(x)W t \^} = J dt (^|W;A(M 3 )W^) = (9a) 
o o 

oo 

J dt(V|iy;(-iA(M 3 ) + {H- |A(M 3 ) - = (9b) 

dt {ip\(w;w t + w;w t )W) = (9c) 





oo 





oo 

dt^\W:WM) = (V«Wo|V> = m) = 1 , (9d) 

provided Wt — > as t — > oo, which is the case if the spectrum of A(M 3 ) is bounded 
away from zero (as in the GRW case). 1 The same calculation for a time integral from 
to t shows that the probability of a flash before t equals 1 — HW^H 2 ; in particular we 
can see why W t should not be unitary. In the GRW case, since A(M 3 ) is a multiple of 

1 Otherwise, the theory is not ill-defined, but rather the flash rate is so low for some ip that there is 
positive probability that no flash occurs. 
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the identity and thus commutes with H, we find that the exponential (jHJ) splits into a 
product of two exponentials, 

W t = e-t/^Q-im/h for t > o . (10) 

The joint distribution of the first n flashes is defined to be 

P(Xi G d 3 ^!, Ti G dt u . . . , X n G d 3 a;„, T n G dt n ) = 

\\K n (0, sex, , aj n , t„) ^|| 2 d 3 «i dti • • • d 3 x n dt n , (11) 

where K n is an operator- valued function on (space-time) n defined by 

A^) 1 / 2 W t „_ t „_ 1 A( a;n _ 1 ) 1 / 2 ^„_ 1 _ 4n _ 2 • • ■ A(a ;i ) 1/2 . (12) 

The square-roots exist since the A (a?) are positive operators. Observing that, by a 
reasoning analogous to Q, 

oo 

J d*x n J dt n K* K n = K*_, K n _ x , (13) 

tn-l 

we see two things: firstly that the right hand side of (jlljl is a probability distribu- 
tion on (space-time)™, and secondly that these distributions, for different values of n, 
are marginals of each other, thus forming a consistent family and arising from a joint 
distribution of infinitely many random variables X i, T\, X 2 , T 2 , . . .. 

The original GRW model contains one further complication: that collapses can act 
on different coordinates, as encoded in the particle index i in (|2*| . Reflecting the fact that 
the model should account for the quantum mechanics of A" distinguishable particles, we 
simply postulate that there are A" different types of flashes, or, equivalently, that each 
flash is labeled by an index i G {1, . . . , N}. Correspondingly, we need to be given A" 
positive-operator- valued functions Ai(x), while Jf? and H are the same for all types of 
flashes. Thus, with every flash (X k ,T k ) is associated a random label I k G {1, . . . , AT}, 
and the joint distribution is defined to be 

F(X 1 G d 3 ^!, T 1 edt 1) I 1 = i 1 ,...,X n e d 3 x n , T n G dt n , I n = i n ) = 

\\K n (0, X\, t\, ii, . . . , x n , t n , i n 

)V|| d^dti-'-d^dtn, (14) 

where K n is now an operator-valued function on [(space-time) x {1, . . . , A^}]™ defined 
by 

M*») 1/a ^-^Vi^n-i) 1 / 2 WVx-*.-, " " " An(^i) 1/2 W*-* , (15) 
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with 

N 



W t = exp(-§ Ai(M 3 )t - for t > , W t = for t < . (16) 

8=1 



In the same way as before, one checks that 

N 



J d 3 *n J dt n K* n K n = K* n ^K n . 

in— 1 / 

un—l 



(17) 



and 

JV 



^ y d 3 ^! y d*i k* k x =i 

H = l n 



implying that ((H)) is a consistent family of probability distributions. This completes our 
definition of a theory from , H, and Ai(x), . . . , Ajy(aj), including the GRW model. 



3 QFT 

Since can be taken to be the Hilbert space of a QFT and H its (regularized) Hamil- 
tonian, we get a collapse version of that QFT as soon as we have the flash rate operators 
A (a;). The model is then defined by (|Sjl. (jllj) . and (J12)) . If the QFT contains several 
species of particles, we may wish to introduce one type of flash for every species, and 
use eq.s (JHJ), (JHJ), and (fTfij) instead. 

So what would be a natural choice of A(x)l Since the operator A(5) determines 
the flash rate in B, it should represent the amount of matter in B, smeared out by 
a Gaussian with width a/y/2. Two natural choices are the particle number density 
operator N(y) and the mass density operator M(y), that is, 

A(x) = J d 3 yAfe- (x - y)2/a2 N(y) , (19) 

respectively 

A(se) = y d 3 y M e - {x ~ y)2/a2 M{y) , (20) 

with M. a suitable constant. The operators ([19)1 have been considered already in J7| for 
CSL in connection with identical particles, albeit apparently in a role more analogous 
to A(x) l l 2 than to A (a?). The choice of (|2~U|) instead of (fTHjl corresponds to the proposal 
of Pearle and Squires 26 to choose the collapse rate proportional to the mass. 

In general, the number operator N(y) can be expressed in terms of suitable annihi- 
lation operators a(y) and creation operators a*(y) in the position representation, acting 
on a suitable Fock space Jif, by 

N{y) = a*{y)a(y), (21) 
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where the product involves, when appropriate, summation over spin indices. In a non- 
relativistic QFT, a(y) is simply the field operator at the location y. (Of course, since 
in nonrelativistic quantum theories usually the Hamiltonian does not contain terms cre- 
ating and annihilating particles, we would have to add such terms artificially to the 
Hamiltonian to obtain a model with non-conserved particle number.) 

For several species of particles corresponding to several quantum fields, we thus 
obtain several rate density operators Aj(a?). 

4 Second Quantization 

I would like to describe another way of constructing flash rate operators A (a;) for QFT. 
It will turn out equivalent to (|19j). It is based on the second quantization algorithm for 
forming a (bosonic or fermionic) Fock space Jif from a one-particle Hilbert space MU\, 
and consists of an algorithm for forming flash rate operators A(ac) acting on Fock space 
Jif from flash rate operators Am (a;) acting on J£m. This algorithm in turn is based 
on two procedures, one concerning direct sums of Hilbert spaces and the other tensor 
products. 

On the direct sum M{ © M2 of two Hilbert spaces, each equipped with a positive- 
operator-valued function Aj(cc), i = 1,2, the natural way of obtaining a positive- 
operator-valued function is 

A(x) = Ai(x)eA 2 (aj). (22) 

On the tensor product space M{ ® J^, it is natural to consider the two functions 

Ai(jc) = Ai(as) <g> I 2 and A 2 (x) = h <g> A 2 (x) , (23) 

defining a theory with two types of flashes. A relevant property of this choice is that if the 
two physical systems corresponding to M{ and M2 do not interact, H = H\<g>l2+Ii®H2, 
and if the initial state vector factorizes, if) = if>i ®ip2, then type-1 and type-2 flashes are 
independent of each other. Indeed, the type-1 flashes are also independent of H 2 , A 2 , 
and ip2, and vice versa. Each of the two systems behaves as if it was alone in the world, 
obeying its own version of the law (jlljl . and that is a reasonable behavior. 

Suppose we do not want two types of flashes, but one. Observe that, by (jU, the 
rate (at time t — 0) for a flash of any type to occur in a set B C M? is (ip\Ai(B)\ip) + 
{if)\A 2 {B)\ip) , the same as the rate of flashes of only one type with rate operators 

A(aj) = Ai(aj) ®I 2 + h ® A 2 (x) , (24) 

which is thus a natural choice of a positive-operator- valued function on the tensor prod- 
uct. 
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On the iV-th tensor power ^® of the one-particle Hilbert space Jf?(i), the flash 
rate operator corresponding to (|24)h 

N 

A(x) = ® A(i)(aj) ® I®^ -0 , (25) 

i=l 

can be written using the permutation operators U a on J#^® N for cr in the permutation 
group S N , 

A ( x ) = m=W l E ^(Ad)^)®/^- 1 ^., (26) 

and therefore assumes values in the symmetric operators, mapping in particular sym- 
metric (bosonic) vectors to symmetric ones and anti-symmetric (fermionic) vectors to 
anti-symmetric ones, thus defining two positive-operator- valued functions A^(aj) act- 
ing on the bosonic (+) respectively fermionic (— ) iV-particle Hilbert space: A^ N Jx) is 
the restriction of the A (a;) given by to the bosonic respectively fermionic subspace 
of jtf* N . The flash theory for N bosons or fermions with these flash rate operators 
is closely related, in a way that will become clear in the next section, to the collapse 
process proposed by Dove and Squires [TTj IT2"]. 

Adding the A^ N Jx) in the sense of (122)) from = to oo yields two functions 
A ± (a;) acting on the bosonic respectively fermionic Fock space. This completes our 
construction for the "second quantization" of Am (a;). If we take Am (a;) to be the 
multiplication operator (J2J) with N = 1 (and % = 1) and a(x) the canonical annihilation 
operator on Fock space, then A.(x) coincides with (|19|) . 



5 Collapses 

After talking so much about flashes, I should point out what they have to do with 
collapses of the wave function. Suppose that n flashes have occurred between time 
and time t, with the fc-th flash at time tk and location x^. Then the distribution of the 
next m flashes after time t, conditional on the history of flashes between and t, is, as 
a consequence of (JTTJ), given by 

IP^-X"n+i £ d 3 a? n+ i, T n+ \ E dt n+ i, . . . , X n+m E d s x n+rn , T n+m E dt ni 
X\ = Xi, T\ = t\, . . . , X n = x n , T n = t n < t, T n+ \ > t 
|| A" m (t, x n _i_i, t n _i_i, . . . , x n + rn , ^n+m)^!! d x n j r \ dt n _i_i ■ d x n j rm dt n _i_ m , (27) 

where ipt is the conditional wave function 

, Wt-tn K n (0, Xx, ti, ■ ■ ■ , x n , t n ) j> 

Pt ~ \\W t . tn K n (0,x 1 ,t 1 ,...,x n ,t n )^W [ ' 
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As a corollary, the flash rate at time t in a set B C M 3 is 

(^|A(S)|^). (29) 
It is ^ that collapses whenever a flash occurs, say at (X,T), according to 

fe + = iiA(x)4 T -ir (30) 

and evolves deterministically between the flashes according to the operators Wt (up 
to normalization). More explicitly, for the example case of N identical particles, the 
(anti)symmetric wave function ip(ri, . . . , r^) collapses to (a normalization factor times) 



/ 

^ i=i 



N_ v 1/2 

e -(x-n)'/«M VCri,...,^) (31) 



for bosons (+) respectively fermions (— ). It is this formula for the collapsed wave 
function that Dove and Squires proposed [T2"] . 

Since the conditional wave function %p t is random, one can form the density matrix 

pt= f P(^Gd0)|0)(0| (32) 



of its distribution, in other words the density matrix of an ensemble of systems, each 
of which started with the same initial wave function ip but experienced flashes inde- 
pendently of the other systems. For the sake of completeness we note that it can be 
computed to be 

„ t t 

Pt = J2j d 3 ^---d 3 x n y dt a ... J At n W t . tn K n \^)^\KlWl tn (33) 
n=0 o 

with K n = K n (0, xx, t\, . . . , x n , t n ), and obeys the master equation 



dt 



i[H,p t ) - i{A(M 3 ), Pi } + J d 3 xA(xy/ 2 p t A(xy/ 2 (34) 



with [ , ] the commutator and { , } the anti-commutator. 

It is tempting to regard the collapsed wave function ip t as the ontology, but I insist 
that the flashes form the ontology. This is a subtle point. After all, since for example 
the wave function of Schrodinger's cat, (| alive) + | dead)) /a/2, quickly collapses into es- 
sentially either | alive) or |dead), it may seem that the collapsed wave function represents 
reality. However, that this view is problematic becomes evident when we note that even 
after the collapse into ipt ~ | alive), the coefficient of |dead) in ip t is tiny but not zero. 
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How small would it have to be to make the cat alive? The more fundamental problem 
with this view is that while the wave function may govern the behavior of matter, it is 
not matter; instead, matter corresponds to variables in space and time j2], called "local 
beables" by Bell jS] and "primitive ontology" by Diirr, Goldstein, and Zanghi [T2j . 

In what I described in the previous sections, the flashes form the primitive ontology. 
But other choices are possible, and this fact underlines that the theory is not completely 
specified by the stochastic evolution law for ip t alone. An example of a different primitive 
ontology, instead of flashes, is the matter density ontology, a continuous distribution of 
matter in space with density 

m(aj,t) = (^|A(aj)|^) (35) 

in our notation. While the two theories (using the same wave function with either 
the flash ontology or the matter density ontology) cannot be distinguished empirically, 
they differ metaphysically and physically. For example, the equations I considered in 
j2H] define a Lorentz-invariant theory with the flash ontology but not with the matter 
density ontology, and strong superselection rules can hold with the flash ontology but 
not with the matter density ontology. For further discussion of the concept of primitive 
ontology see 

6 Predictions 

The empirically testable predictions of the collapse QFT model we described agree with 
the standard predictive rules of QFT to the same extent as the GRW model (say, with 
mass-dependent collapse rate) agrees with standard quantum mechanics. To see this, 
recall first from the paragraph containing eq. (J23*j) that when a system (defined, e.g., by 
a region in 3-space) can be regarded as decoupled and disentangled from its environment 
then its flash process is independent of the environment. Next note that, by (|29|). the 
total flash rate is (^ t \A(M. 3 )\^ t ) , proportional to either the (value regarded in quantum 
theory as the) average number of particles (relative to ipt) or the average net mass. 
As a consequence, as with the GRW model, a system containing fewer than a thousand 
particles experiences no more flashes than once in 100,000 years. Up to the first flash, the 
deviation of tp t from the Schrodinger evolution is small for t <C t/ AN if AN is a bound 
on the spread in particle number of e~ lHs ^ h ip, < s < t. A macroscopic piece of matter, 
in contrast, with over 10 22 particles, experiences millions of flashes every second. A 
macroscopic superposition essentially breaks down with the first flash (with consequences 
for the distribution of the future flashes) to one of the contributions, and for the same 
reasons as in the GRW model the random choice of the surviving contribution occurs 
with almost exactly the quantum theoretical probabilities. And like in the GRW model, 
this entails that experiments that "measure" any quantum "observable" on a microscopic 



10 



system have almost exactly the quantum theoretical distribution of outcomes. Dove 
and Squires [12] provide a discussion of the consistency of their collapse process for 
identical particles with those experimental data that yield restrictions on the possibility 
of spontaneous collapse. 

Some collapse theories imply the possibility of superluminal (i.e., faster than light) 
signalling; even if the theory is hard to distinguish empirically from standard quantum 
theories, those experiments sensitive enough to detect the deviation can allow signalling 
using EPR-Bell pairs. Such collapse theories are therefore unlikely to possess a Lorentz- 
invariant version. In contrast, the collapse QFT developed here and the GRW model 
exclude superluminal signalling; this follows essentially from their property that the 
distribution of the flashes is quadratic in if). Indeed, if two systems are entangled but 
decoupled, H — Hi g) I 2 + I\ ® H 2 , and the flash rate operators are additive according 
to or (jUJ), then, as a consequence of (fTTj) . the marginal distribution of the flashes 
belonging to system 1 depends on H, A, and if) only through Hi, Ai, and the reduced 
density matrix tr 2 \if>) {if>\ of system 1. To see this, consider first the case (|23|). in which 
there are two types of flashes for the two systems; in this situation the claim follows from 
the fact that the operators K n defined by (fT5|) decompose into Ki m ®K 2)1l2 . In the case 
(J2*3jl of a single type of flashes it is not obvious which flashes are to be attributed to which 
system, unless the two systems have disjoint supports (Si,S 2 C M 3 with Si H 52 = 0) 
and we count the flashes in Si for system 1. But then we can in fact regard the flashes as 
labeled, the label being a function of the location, corresponding to Aj(jc) = lsi( x ) A{ x ) 
with lg. the characteristic function of Si, which brings us back to the previous case. 

7 Literature 

GianCarlo Ghirardi sometimes suggests in his writings that identical particles or QFT 
cannot be treated in the framework of the GRW model in a satisfactory way fSJ 
page 118], [3, pages 312 and 382], but require a diffusion process in Hilbert space; I 
think that the model I have presented (respectively, as far as A" identical particles are 
concerned, the model of Dove and Squires ^2]) is a counterexample. Part of the reason 
why the model of Dove and Squires has not received enough attention may be that they 
have not made clear enough, in my view, its naturalness and simplicity, and that they 
have presented it on equal footing with another, much less natural, proposal. 

I know of five variants of the GRW model for identical particles that have been 
proposed, apart from the one discussed in this paper: One was introduced by Ghirardi, 
Nicrosini, Rimini, and Weber in 1988 ^H], which, however, appears theoretically unsat- 
isfactory since it prescribes that the flashes of a system of A" identical particles occur 
in clusters of A" simultaneous ones, leaving no hope for a Lorentz-invariant version. It 
is also presumably empirically inadequate since it predicts that a superposition of two 
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wave packets for the same single particle at a distance greater than 10 kilometers col- 
lapses within 10~ 7 sec. The second proposal was made by Kent in 1989 (20]; i n which, 
however, the distribution of the flashes is not quadratic in tp, thus allowing superluminal 
signalling. The same is true of the third proposal, which is contained as well in the pa- 
per of Dove and Squires ^2]- The fourth proposal was made by Bell in 1987, but never 
published, except in a brief description in Sec. IV C of ^7]. I have learnt about it from 
Alberto Rimini on a recent conference in honor of GianCarlo Ghirardi's 70th birthday 
at Trieste and Mali Losinj; there I have learnt as well that equations similar to (jHJ) and 
(J3*U|) had already been considered in a different context, namely for models of coupling 
classical and quantum systems [7j. The fifth variant was mentioned by Ghirardi, Pearle, 
and Rimini in 1990 as a side remark in Sec. IV C of ^7j: The process for the wave 
function is similar in spirit to the proposal of ^H]; involving a simultaneous collapse 
for all particles, but differs in that it cannot be associated with flashes, as the collapse 
involves a continuous distribution function n(x) on, rather than a set of N points in, 
3-space. 

Acknowledgments. I thank Philip Pearle for his detailed comments on a previous version 
of this article. 
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